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Abstract 

The differential cross section of the process 77 7r°7r° has been measured in the kinematical 
range 0.6 GeV <W < 4.0 GeV and | cos 0*\ < 0.8 in energy and pion scattering angle, respectively, 
in the 77 ccnter-of-mass system. We find at least four resonant structures including a peak from 
/o(980). In addition, there is evidence for XcO production. We also make a preliminary discussion 
of the angular dependence and cross section ratio to 77 — > tt+tt". 

PACS numbers: 13.20.Gd, 13.60.Le, 13.66.Bc, 14.40.Cs,14.40.Gx 
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I. INTRODUCTION 



Measurements of exclusive hadronic final states in two-photon collisions provide valuable 
information concerning the physics of light and heavy-quark resonances, perturbative and 
non-perturbative QCD and hadron-production mechanisms. So far, we have measured the 
production cross sections of charged-pion pairs [H, [21, charged and neutral-kaon pairs jl, 0], 
and proton-antiproton pairs We have also analyzed D-meson-pair production to search 
for a new charmonium state [5)]- 

In this report, we show an analysis of neutral-pion pair production in two-photon pro- 
cesses. The motivation for this study is essentially the same as that for the charged-pion 
pair case. But the two processes are physically different and independent; we cannot predict 
very precisely what happens in one by measuring only the other. 

In the low energy region {W < 1.0 GeV), it is expected that the difference of meson 
electric charges plays an essential role in the difference between the tt+tt^ and Tr^Tr*^ cross 
sections. The predictions are not easy because of non-perturbative effects. In the interme- 
diate energy range (1.0 < W < 2.4 GeV), the formation of meson resonances decaying to 
TTTT is the dominant contribution. For ordinary qq mesons conserving isospin in decays to 
TTTT, the only allowed J^*-^ states produced by two photons are 0"'"(even)+"'', that is, /j=even 
mesons. The ratio of the /-meson's branching fractions, B{f TT^TT^)/B{f tt+tt^) is 1/2 
from isospin invariance. But, interference of the resonance with the continuum component 
which cannot be precisely calculated distorts the ratio even near the resonant peaks. The 
TT^Tc^ channel has an advantage in the study of resonances, since a smaller contribution from 
the continuum is expected in it than in the tt+tt" channel. 

For higher energies, we can invoke a quark model. In the leading order calculations [g], 0] 
which take into account spin correlation between quarks, the tt^tt" cross section is predicted 
to be much smaller than that of tt+tt^, and the cross section ratio of vr^vr^ to tt+tt" is around 
0.03-0.06. However, higher-order or non-perturbative QCD effects can modify the ratio. 
For example, the handbag model which considers soft hadron exchange predicts the same 
amplitude for the two processes, and this ratio becomes 0.5 jl]. Analyses of energy and 
angular distributions of the cross sections are essential for determining properties of the 
observed resonances and for testing the validity of QCD models. 

We present preliminary results of the measurement of the differential cross section, 
da/d\cos6*\, for the process 77 — > vr^vr^ in a wide two-photon center-of-mass (cm.) en- 
ergy (W) range from 0.6 to 4.0 GeV, and a cm. angular range, | cos6'*| < 0.8. (Although 
we do not need to put the absolute- value symbol for cos 6* for the present channel where the 
identical particle pairs appear in both initial and final states, we follow the usual conven- 
tion of the two-photon differential cross section to avoid unnecessary confusion.) Our data 
sample is by several hundred times larger than in previous experiments [lo| . 



II. EXPERIMENTAL APPARATUS AND TRIGGER 

We use a data sample that corresponds to an integrated luminosity of 95 fb^^ recorded 
with the Belle detector at the KEKB asymmetric-energy e~^e~ collider [llj. The energy of 
the accelerator was set at 10.58 GeV (83 fb'^), 10.52 GeV(9 fb'^), 10.36 GeV(T(35) runs, 
2.9 fb~^) and 10.30 GeV(0.3 fb~^), in the cm. energy of the e~^e~ beams. The difference of 
the two-photon flux (luminosity function) in the measured W regions due to the different 
beam energies is very small (a few percent at maximum), and the fraction of integrated 



4 



luminosity of the runs with the lower beam energies is very small. We therefore combine the 
results for different beam energies. The effect on the cross section deviation is less than 0.5%. 
The analysis is made in the "zero-tag" mode, where neither the recoil electron nor positron 
is detected. We restrict the virtuality of the incident photons to be small by imposing strict 
transverse-momentum balance with respect to the beam axis for the final-state hadronic 
system. 



A comprehensive description of the Belle detector is given elsewhere jl2| . We mention here 



only those detector components which are essential for the present measurement. Charged 
tracks are reconstructed from hit information in a central drift chamber (CDC) located in 
a uniform 1.5 T solenoidal magnetic field. The detector solenoid is along the z axis which 
points in the direction opposite to that of the positron beam. The CDC and a silicon 
vertex detector measure the longitudinal and transverse momentum components (along the 
z axis and in the r^? plane, respectively). Photon detection and energy measurements are 
performed with a CsI(Tl) electromagnetic calorimeter (ECL). 

In the present measurement, we require that there be no reconstructed tracks coming from 
the vicinity of the nominal collision point. Therefore, the CDC is used for vetoing events 
with track(s). Photons from decays of two neutral pions are detected and their momentum 
vectors are measured by the ECL. 

Signals from the ECL are used to trigger the signal events. The conditions of the ECL 
trigger are as follows: The ECL total energy deposit in the triggerable acceptance region 
(see the next subsection) is greater than 1.15 GeV (the 'HiE' trigger), or, the number of 
the ECL clusters counted according to the energy threshold at 110 MeV for segments of the 
ECL is four or greater (the 'Clst4' trigger). The above energy thresholds are determined 
from the experimental data from a study of the correlations between the two triggers. No 
software filterings are applied for triggering events by either or both of the two ECL triggers. 



III. EVENT SELECTION 

The selection conditions for signal candidates of 77 -n^n'^ are as follows: All the 
variables in the criteria (l)-(5) are measured in the laboratory frame; (1) there is no good 
track that satisfies cir < 5 cm, \dz\ < 5 cm and pt > 0.1 GeV/c, where dr and dz are the 
radial and axial distances, respectively, of the closest approach (as seen in the np plane) 
to the nominal collision point, and the pt is the transverse momentum measured in the 
laboratory frame with respect to the z axis; (2) the events are triggered by the HiE or Clst4 
triggers; (3) there are two or more photons whose energy is greater than 100 MeV; (4) there 
are just two vr^'s, each having transverse momentum greater than 0.15 GeV/c, and each of 
the decay-product photons has energy greater than 70 MeV, where the ir^ is reconstructed 
from two photons and is selected with a value of a mass constraint fit; (5) the total energy 
deposit in ECL is smaller than 5.7 GeV. 

Then, the transverse momentum in the e+e^ cm. frame (|SpJ'|) of the two-pion system 
is calculated. For further analysis, (6) we use events with |Sp^| < 50 MeV/c as the signal 
candidates. 

In order to reduce uncertainty from the efficiency of the hardware ECL triggers, we set 
offline selection criteria which emulate the hardware trigger conditions as follows: (7) the 
ECL energy sum within the triggerable region is greater than 1.25 GeV, or all the four 
photons composing the two 7r° are contained in the triggerable acceptance region. Here, we 
define the triggerable acceptance region as the polar-angle range in the laboratory system 
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17.0° <e < 128.7°. 



IV. SIGNAL CANDIDATES AND BACKGROUNDS 

A. Yield distribution of the signal candidates 

We derive the cm. energy W of the two-photon coUision from the invariant mass of the 
two neutral pion system. We calculate cosine of the scattering angle of 7r° in the 77 cm. 

frame, |cos^^*| for each event. We use the e~^e~ collision axis in the e~^e~ cm. frame as 
the reference of the polar angle as an approximation, because we do not know the exact 77 
collision axis. 

The two-dimensional yield distribution of the selected events is shown by the lego plot in 
Fig. 1. The W distribution with | cos^*| < 0.8 is shown in Fig. 2(a). We find at least four 
resonant structures: clear peaks from /o(980) near 0.98 GeV and /2(1270) near 1.25 GeV, 
and broad structures near 1.65 GeV and near 1.95 GeV. Figure 2(b) is an enlarged view of the 
charmonium region for | cos^*| < 0.4, where we see some hints of charmonium contributions 
in the Xco(~ 3.40 GeV) and Xc2(~ 3.55 GeV) mass region. 

B. Background subtraction 

We study the pt-balance distribution, i.e., the event distribution in |Spj | to separate the 
signal and background components. The signal Monte Carlo (MC) events show that the 
signal component peaks around 10-20 MeV/c in this distribution. In the experimental data, 
however, in addition to such a signal component, we find some contribution from the pt- 
unbalanced components in the low- 14^ region. The source of such p^-unbalanced components, 
in general, is considered to be non-exclusive processes such as 7r°7r°7r° etc. But, in the present 
case, the background found in the experimental data is very large only in the low W region 
where the vr^'Tr^'Tr^ contribution is expected to be much smaller than 7r°7r° in two-photon 
collisions. (Note that a C-parity odd system cannot go to tt^tt".) We believe that the 
backgrounds are dominated by beam-background photons (or neutral pions from secondary 
interactions) or spurious hits in the detector. 

Figures 3 (a) and (b) show the p^-balance distributions in the low W region. With 
the fit described below, we separate the signal components from the background. In the 
intermediate or higher energy regions, the backgrounds are less than 1%, buried under the 
/2(1270) peak (Fig. 3(c)), or we do not see any detectable pt-unbalanced components within 
the statistical errors (as shown in Fig. 3(d)). Even at the highest energy 3.6 GeV< W < 
4.0 GeV, I cos^*| < 0.4, we find no visible contamination from the background. 

A fit to the prbalance distribution is made to separate the signal and background com- 
ponents for the W region below 1.2 GeV. The fit function is a sum of the signal component 
and the background component. The signal component is taken as an empirical function 
reproducing the shape of the signal MC, y = Ax/i^x"^'^ + B + Cx), where {x = |SpJ'|, A, B 
and C arc the fitting parameters, and y is the distribution). This function has a peak at 
X = (j-[)~ and vanishes at x = and infinite x. The shape of the background is taken as 
a linear function y = ax for x < 0.05 GeV/c which is smoothly connected to a quadratic 
function above x > 0.05 GeV/c. 
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FIG. 1: The two-dimensional distribution of the experimental candidates. The same distri- 
bution is viewed from two different directions. 

The background yields obtained from the fits are fitted to a smooth two-dimensional 
function of (W, \ cos^*|) in order to minimize statistical fluctuations. Then, the backgrounds 
are subtracted from the experimental yield distribution. The background yields integrated 
over the angle is shown in Fig. 2(a). We omit the data points in the small-angle (| cos^*| > 
0.6) region in W < 0.72 GeV, because there the background dominates the yield. 



C. Unfolding the W distributions 

We estimate the invariant-mass resolutions from studies of the signal-MC and experi- 
mental events. We find that the MC events show a relative invariant-mass resolution of 
1.4%, which is almost constant in the whole W region of the present measurement. The 
momentum resolution of 7r° is known to be worse by about 15% in the experimental data 
than in MC from a study of the p^-balance distributions (described in Sect. V). Moreover, 
the distribution of the MC is asymmetric; it has a longer tail on the lower mass side. 

An asymmetric Gaussian function with standard deviations of 1.9%!^ and 1.3%iy on 
the lower and higher sides of the peak, respectively, is used and approximates the smearing 
reasonably well. 

We calibrate the experimental energy scale and invariant-mass distribution using the 77 
invariant mass from experimental samples of rj' —>■ 77 from two-photon processes. The peak 
position is consistent with the nominal mass of rj' with an accuracy better than 0.2%. 

This invariant- mass resolution is comparable to or larger than the W bin width (20 MeV) 
used in Figs. (1) and (2). We made an unfolding of the invariant-mass distribution in each 
I cos 6* I bin separately, to correct for migrations of signal yields between different W bins, 
based on the smearing function of the above asymmetric Gaussian shape. The migration in 
the I cos 6** I direction is expected to be small and is neglected. 

The unfolding is made using the Singular Value Decomposition (SVD) algorithm 13|] in 
the yield level, and is applied so as to reproduce the corrected W distribution in the 0.9 - 
2.4 GeV region, using data at observed W between 0.72 and 3.0 GeV. For lower energies, 
W < 0.9 GeV, the effect of the migration is expected to be small because of the smaller 
invariant-mass resolution compared with the bin width. For higher energies, W > 2A GeV, 
where the statistics is relatively low and the unfolding would enlarge the errors, we rebin 
the result with a bin width of 100 MeV, instead of unfolding, as described in Sect. IV.C. 
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FIG. 2: The W distribution of the candidate events, (a) | cos6*\ < 0.8. The curve is an estimate 
of backgrounds from events with pt imbalance, (b) |cos0*| < 0.4, near the charmonium region. 
The curve is the fit described in Sect. VIII. 

Distributions before and after the unfolding for a typical angular bin (| cos^*| = 0.225) are 
shown in Fig. 4. 

V. DETERMINATION OF EFFICIENCY 

We determine the trigger efficiency for signal events using the detector and trigger simu- 
lators applied to the signal MC events. ^ 

The signal MC events for e^e^ —>■ e^e^n^Ti^ are generated using TREPS code [14]] for 
the trigger efficiency study at 27 fixed W points between 0.5 and 4.1 GeV, isotropic in 
I cos 6* I . The angular distribution at the generator level does not play a role for the efficiency 
determination, because we calculate the efficiencies separately in each | cos 6* \ bin with width 
0.05. 5 X lO'' events are generated at each W point. 
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|Sp^*| (GeV/c) 



FIG. 3: The distribution of imbalance in pt for candidate events, (a) In the bin centered at 
W = 0.90 GeV and | cos 6**1 = 0.05 (The bin width is 0.04 GeV and 0.1 in the W and jcos^l 
directions, respectively, in (a)-(c).), the experimental distribution (histogram) is fitted with the sum 
of signal and background components (curves). The grey region shows the estimated background 
contamination in the signal region, (b) W = 0.66 GeV, same as in (a) for the others, (c) In the bin 
oiW = 1.18 GeV, | cos^*| = 0.65, the experimental distribution (dots with error bars) is compared 
with the signal MC (histogram). The grey region shows the estimated background contaminating 
the signal region obtained from the fit. (d) For 3.6 < W < 4.0 GeV and |cos0*| < 0.4, the 
experimental distribution (dots with error bars) is compared with the signal MC (histogram). 

To minimize statistical fluctuation in the MC, we fit the numerical results of the trigger 
efficiency to a two-dimensional empirical function in {W, \ cos^*|). The W dependences are 
shown in Fig. 5 for two typical angular bins. We find that the trigger efficiency is almost 
fiat and close to 100% for the region W above 1.3 GeV. The angular dependence is rather 
small. 

Separately, we generated signal MC events at 48 W points in the same W region with 10^ 
events at each value of W, for the acceptance calculations. Here we call the efficiency of the 
selections not including the hardware trigger "acceptance" ; the net efficiency is a product of 
the trigger efficiency and the acceptance. The determined acceptance from the MC events 
is fitted by a smooth two-dimensional function of (W, \ cos 6**1) (Fig. 6). It is about 11% at 
maximum and gets smaller (down to around 1%) at lower W or smaller cm. angle (larger 

|cosr|). 

The acceptance calculated using the signal MC events is corrected for a systematic dif- 
ference found between the peak widths in the p^-balance distributions of the experimental 
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FIG. 4: Invariant mass distributions of the yield in each bin before(orange colored triangles) and 
after (dark-blue diamonds) the unfolding, at |cos^*| = 0.225. The inset shows a logarithmic plot 
of the same distribution. The bin width is different for the regions W below and above 2.0 GeV. 




FIG. 5: The W dependence of the trigger efficiency for two angular bins. The curves are the fit to 
parameterize it. 



data and the MC, which could affect the acceptance through the |Sp*| cut. It originates 
from a difference in the momentum resolution for vr" between data and MC. We find that 
the peak position of the experimental data in the distribution is 10% to 20% higher than 
the MC expectation, depending on W and | cos6'*|. The correction factor ranges from 0.90 
to 0.95. 
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FIG. 6: The W dependence of the acceptance for two angular bins. The curves are the fit to 
parameterize it. No correction for the vr*^ momentum resolution is included. 

VI. CROSS SECTION CALCULATION 

The differential cross section for each {W, |cos6'*|) point is derived by the following 
formula: 

da _ AY - AB 

d\ cos9*\ ~ AWA\ cos9*\ J CdtL^^{W)r]trg7]^cc 

where AY and AB are the signal yield and the estimated p^-unbalanced background in the 
bin, AW and A| cos 6**1 are the bin widths, / Cdt and L^^{W) are the integrated luminosity 
and two-photon luminosity function calculated by TREPS, respectively, and r^trg and ?7acc 
are the trigger efficiency and the acceptance, respectively, the latter including the correction 
described in the previous section. 

Figures 7(a) and (b) show the W dependence of the cross section integrated over | cos 9* \ < 
0.8 and | cos 6**1 < 0.6, respectively. They are obtained by simply adding da / d\cos9*\ ■ 
A I cos^*| over the corresponding angular bins. 

The data points for 0.9 GeV < W < 2.4 GeV are the unfolded results with the bin 
widths AW = 0.02 GeV (0.04 GeV) for W above (below) 2.0 GeV. For the data points 
above 2.4 GeV, we average five data points each with a bin width of 0.02 GeV and get 
results for every 0.1 GeV at W. We have removed the bins in the range 3.3 GeV < W < 
3.6 GeV, because we cannot separate the components of the Xco, Xc2 and the continuum in 
a model-independent way due to a finite mass resolution and insufficient statistics in the 
measurement. 

We show the angular dependence of the differential cross section at several W points in 
Fig. 8. 

It is noted that the cross section results after the unfolding are no longer independent of 
each other in the neighboring bins, in both central values and size of errors. 
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FIG. 7: The cross section results integrated in the angular regions (a) |cos^*| < 0.8 and (b) 
I cos 0*1 < 0.6 after the unfolding and rebinning. 



W - 0.97 GeV 



W= 1.27 GeV 



50 






AS 






AO 






35 






30 






25 




20 






15 






to 






5 














W = 1.95 GeV 



0.2 OA 0.6 0.8 

Icose-i 




0.2 0.4 0.6 0.6 
IC08 9-I 



250 
200 
150 
100 
50 

oL 



0.2 0.4 0.6 0.8 1 

W= 3.15 GeV 




0.2 0.4 0.6 0.8 

Icose*! 




0.2 0.4 0.6 0.6 
|a»e*| 



FIG. 8: The differential cross sections for five selected W points, 0.97 GeV, 1.27 GeV, 1.95 GeV, 
2.45 GeV and 3.15 GeV. The results of the first three W points are after the unfolding. 



VII. SYSTEMATIC ERRORS 

We summarize the evaluation of the systematic errors for a{\ cos^*| < 0.8) {a{\ cos6*\ < 
0.6) with W < 0.72 GeV) for each W point. They come from the following major error 
sources: 

Trigger efficiency: The systematic error of the Clst4 trigger is assigned as 2/3 of the dif- 
ference of the efficiencies with the different threshold assumptions for the ECL cluster - 
110 MeV and 100 MeV - set in the trigger simulator for the energy region IV < 2.5 GeV. 
Separately, we take the uncertainty in the HiE trigger efficiency as 4% for the whole W re- 
gion. The systematic errors from the two triggers are added in quadrature. This systematic 
error is large in the low W region, 20%-30% for W < 0.8 GeV. 
The reconstruction efficiency: 6% for two pions. 

The pt-balance cut: 3%-5%. One half of the correction size discussed in Sect. V. 
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FIG. 9: The cross section results integrated over the angular regions | cos0*| < 0.8 compared with 
the previous measurement from the Crystal Ball [lo| . 



Background subtraction: 20% of the size of the subtracted component is assigned as the 
error from this source. In the W region where the background subtraction is not applied 
{W > 1.2 GeV), we neglect the error for 1.2 GeV< W < 1.5 GeV, or assign 3% for 
ly > 1.5 GeV, considering an upper limit of the background contamination expected from 
the pt-balance distributions. 

Luminosity function 4-5%. 4% (5%) for W < (>)3.0 GeV. 

The total systematic error is 9% in the intermediate W region, 1.05 GeV< W < 2.7 GeV. 
It is much larger for lower W, 15% at W = 0.85 GeV, 30% at W = 0.70 GeV and 55% 
at W = 0.61 GeV. The error is dominated by the background subtraction for low W. For 
higher W, the systematic error is rather stable, 10%-11% for 2.7 GeV< W < 4.0 GeV. 



VIII. DISCUSSION 

We compare our results with the previous measurements by Crystal Ball at DORIS 



II [lOj (Fig. 9). The agreement is fairly good. The error bars from the two experiments are 
statistical only, and the systematic errors ( 7% (11%) for W > (<) 0.8 GeV for the Crystal 
Ball results) should also be considered in the comparison. The present measurement has 
several hundred times more statistics than the Crystal Ball measurement. 

We find prominent resonant structures near 0.98 GeV and 1.27 GeV. They are from 
/o(980) and /2(1270), respectively and are observed also in the 77 tt+tt" process [l[. It 
is for the first time that the /o(980) is observed as a clear peak in 77 ti^tt^. We find clear 
evidence for rather broad peaks at 1.65 GeV and 1.95 GeV. For these, any quick assignment 
to well known states is not easy. In addition, we find a hint of a possible inflection point 
near 1.4 GeV and a structure in the mass region of two charmonium states Xco and Xc2- The 
structures above found in the 1.2-2.1 GeV region are somewhat similar to the distribution 



observed in the vr'^Tr^ spectrum from the 7r~p — >■ n^n^n experiment, GAMS 15 . 

We fit the yield distribution in the range 2.8 GeV< W < 4.0 GeV and | cos 6^*1 < 0.4 in- 
cluding the contribution from the Xco and Xc2 charmonia with a binned maximum-likelihood 
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FIG. 10: The cross sections of the 77 — > it tt and 77 — > 7r"''7r~ reactions for | cos^*| < 0.6. The 
blue closed circle and the violet triangles are the cross sections of the two reactions in units of nb. 
The red closed circle is the ratio. The error bars are statistical only. 



method. The fit is shown in Fig. 2(b). We take into account the finite invariant-mass reso- 
lution effect introduced in the last section in the fit. We have fixed the nominal masses and 
the width of Xco to the world averages jl6| (we neglect the width of Xc2)- The background 
component is assumed to have the shape of ~ VT"". 

From the fit, we find the yields 35.3 ± 9.2 events and 8.2 ± 6.4 events for Xco and 
Xc2, respectively. These yields provide the products of the two-photon decay widths and 
the branching fractions, Tyy{xcj)B{Xcj — 7r°7r°) = 8.4 ± 2.2{stat.) ± 0.8{syst.) eV and 
0.29 ± 0.23 ± 0.03 eV, for Xco and Xc2, respectively. The former and latter provide ev- 
idence for the 77 — s> Xco ~^ ti^tt^ signal at the 4.5(T level and an upper limit at the 
90% C.L., T^^{xc2)B{Xc2 7r°7r°) < 0.75 eV which is obtained from the yield where 
the two times the log-likelihood of the fit is smaller by (1.64)^ than that of the best fit, 
respectively. The above central values show good agreement with the world averages of 
r^^{Xcj)B measured in the Xcj —>■ vr+vr^ process [l6j, considering isospin invariance (apply 
n^Ti^ : TTTT = 1 : 3 to the branching fractions), which includes the Belle measurements of the 
process 77 Xcj tt+tt" [;2j]. 

The general trend of the angular dependence of the differential cross section is as follows: 
The differential cross section has a maximum at | cos6'*| = for the W < 2.1 GeV region. 
For the ly > 1.9 GeV, however, the angular dependence shows a rise toward the forward 
angles. In the higher W region, the point in cm. angle at which the rise in the differential 
cross section begins moves toward the forward direction, and the rise gets steeper, as W 
increases. 

We show the cross section ratio between 77 tt^tt^ and 77 — * tt+tt" j2| for | cos^^*| < 0.6 
and 2.4< W < 4.0 GeV, in Fig. 10. The error bars in the figure are statistical only. Each 
measurement of the cross section has typically a 10% systematic error, partially correlated 
with each other for different points. Several 77 vr+vr" measurements above the charmo- 
nium masses have larger systematic errors, ~ 25%. 

The energy dependence of the ratio above W > 2.7 GeV is smooth with a value around 
0.3-0.4. This ratio is larger than the prediction from the lowest-order QCD calculation. We, 
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however, need more detailed investigations of the W and angular dependences in comparison 
with the predictions of the QCD models to test the expected asymptotic nature. 



IX. CONCLUSION 

We have measured the cross section of the process 77 —>■ 7r°7r° in the 77 cm. energy 
and angular regions of 0.60 GeV < 1^ < 4.0 GeV and | cos^*| < 0.8. In the cross section, 

several resonant structures are seen, including a statistically significant peak from the 
/o(980). We find that the angular dependence (forward- and/or large-angle enhancements) 
changes drastically at around W = 2.0 GeV. The ratio of cross sections of 77 — > 7r°7r° and 
77 7r''"7r~ in the 3 GeV region is also obtained. 
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